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ABSTRACT 


Two-stage  Gijford-McMahon  refrigerators  are  candidates  for  use  in  cooling  superconducting 
magnets  for  naval  applications  in  mine  countermeasures  and  electric-drive  propulsion  for  ships.  For 
these  applications,  the  refrigerators  can  be  expected  to  undergo  shock  and  vibration  due  to  the  mo¬ 
tion  of  the  platform  on  which  they  are  mounted  and  to  explosions  of  nearby  mines.  If  the  refrigerator 
is  to  continue  operating  effectively,  the  cylinder  walls  must  not  be  permanently  deformed  when  sub¬ 
jected  to  shock  loads,  and  therefore  stresses  in  the  walls  must  not  approach  the  elastic  limit  of  the 
wall  materioL  The  stress  in  the  cylinder  walls  due  to  specified  shocks  was  determined  by  an  axisym- 
metric  finite  element  shock  analyses  of  the  two-stage  cylinder  and  displacer  assembly.  For  this  analy¬ 
sis,  it  was  assumed  that  the  displacers  were  at  the  bottom  of  their  stroke  and  each  was  resting  on  the 
bottom  of  its  cylinder.  Constant  horizontal  and  vertical  accelerations  oflOOg  and  a  time-dependent 
acceleration  with  maximum  amplitude  of  103 g  were  applied  to  the  model  The  analysis  for  the  verti¬ 
cal  shock  loading  produced  a  maximum  stress  of  36. 7  MPa,  5  percent  of  yield  for  the  304-type  stain¬ 
less  steel  used  for  the  cylinder  walls.  The  horizontal  shock  loading  produced  a  maximum  stress  of 
144  MPa,  19  percent  of  yield  for  the  stainless  steel  The  time-dependent  analysis  produced  a  maxi¬ 
mum  stress  of  18.4  MPa,  2  percent  of  yield  for  the  stainless  steel  These  stress  values  are  consistent 
with  the  expectation  that  a  constant  acceleration  will  produce  larger  stresses  than  a  time-varying  ac¬ 
celeration  of  equal  maximum  amplitude  that  does  not  excite  a  resonance  of  the  structure.  Given  the 
anticipated  values  of  horizontal  and  vertical  shock  to  be  applied  during  testing,  it  is  expected  that  the 
refrigerator  will  not  experience  damaging  strain. 

ADMINISTRATIVE  INFORMATION 

This  work  was  performed  as  part  of  the  Advanced  Lightweight  Influence  Sweep  System  Advanced  Technology 
Demonstration  Program.  The  program  manager  is  the  Program  Executive  Office,  Mine  Warfare  (PMO-407).  The 
work  was  an  element  of  the  Refrigerator  task  of  the  ALISS  Magnetic  Subsystem  Program  in  the  Electrical  Machinery 
Technology  Branch,  Machinery  Research  and  Development  Directorate.  The  work  was  performed  under  work  unit 
number  1-8120-400  by  the  Computational  Signatmes  and  Structures  Group,  Ship  Systems  and  Programs  Directorate 
(Code  2042). 


INTRODUCTION 

The  U.S.  Navy  is  investigating  applications  of  superconducting  magnet  technology  for  use  in  mine  countermea¬ 
sures  and  electric  drive  propulsion  for  ships.  Two-stage  Gifford-McMahon  refrigerators  are  likely  candidates  for 
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Figure  1.  Balzers  Model  UCH-130  Gifford-McMahon  refrigerator 
cooling  the  magnets  in  these  applications.^  *  Figure  1.  shows  a  photograph  of  a  Balzers  Model  UCH-130  Gifford- 
McMahon  refrigerator  of  the  type  analyzed  here.  For  each  stage  in  this  Gifford-McMahon  design,  the  regenerative  heat 
exchanger,  subsequently  referred  to  as  a  regenerator,  is  located  inside  the  displacer  resulting  in  the  addition  of  non- 
structural  mass  to  the  displacer. 

For  naval  applications,  these  refrigerators  can  be  expected  to  undergo  shock  and  vibration  due  to  the  motion 
of  the  platform  on  which  they  are  mounted.  In  mine-sweeping  applications,  nearby  mine  explosions  may  produce 
additional  shocks.  If  the  refrigerator  is  to  continue  operating  effectively  after  experiencing  shock  loads,  the  cylinder 
walls  must  not  be  permanently  deformed,  and  therefore  the  wall  material  must  not  approach  its  elastic  limit  when 
subjected  to  these  loads.  The  goal  of  this  work  is  to  determine  the  stress  in  the  cylinder  walls  due  to  expected  shocks 
by  making  shock  analyses  of  the  cylinders  of  the  two-stage  cylinder  and  displacer  assembly.  The  cylinders  and  displacers 
were  modeled  with  finite  elements  and  constant  100  g  horizontal  and  vertical  accelerations  and  a  time-dependent 
acceleration  having  a  103  g  maximum  amplitude  were  applied  to  the  model.  The  resulting  displacements  and  stresses 
were  plotted,  and  the  maximum  values  of  each  were  extracted. 

This  report  documents  the  modeling  assumptions,  details,  and  results  produced  by  the  analyses.  Appendix  A 
contains  diagrams  of  the  refrigerator  components  represented  by  the  finite  element  model  and  Appendix  B  contains 
representative  displacement  and  stress  plots. 

*  The  references  are  listed  on  page  23. 
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MODELING  ASSUMPTIONS 


A  two-stage  Gifford-McMahon  refrigerator  mounted  in  a  vacuum  vessel  was  analyzed  using  the  finite  element 
method  to  determine  its  mechanical  response  to  shock  loading.  Since  determining  the  effect  of  shock  on  the  vacuum 
vessel  is  outside  the  scope  of  this  work,  the  prescribed  shock  acceleration  was  assumed  to  be  that  transferred  to  the 
refrigerator  by  the  vessel  and  was  applied  at  the  top  of  the  first  stage  where  the  refrigerator  is  attached  to  the  vessel. 
To  obtain  the  effect  of  a  constant  acceleration  produced  by  a  force  applied  at  the  top  of  the  first  stage,  the  motion  of 
the  top  of  the  first  stage  was  constrained  and  a  gravitational  acceleration  in  the  direction  opposite  to  the  shock 
acceleration  was  applied  to  the  structure.  Figure  3.  in  Appendix  A  shows  the  location  at  which  the  acceleration  loading 
was  applied.  Displacements  and  stresses  produced  by  vertical  accelerations  parallel  to  the  geometric  axis  of  the  cylinders 
and  horizontal  accelerations  perpendicular  to  the  axis  were  determined  using  finite  element  analyses.  The  vertical  and 
horizontal  accelerations  were  applied  in  separate  computer  runs.  Since  the  assembly  of  the  first-  and  second-stage 
cylinders  and  displacers  of  the  refrigerator  is  rotationally  symmetric  with  respect  to  the  central  axis,  the  assembly  was 
modeled  using  axisymmetric  finite  elements.  A  three-dimensional  axisymmetric  structure  is  generated  by  rotating  a 
two-dimensional  region  around  an  axis  of  symmetry  and  an  axisymmetric  analysis  models  the  structure  by  subdividing 
this  two-dimensional  region.  Since  the  subdivided  region  was  two-dimensional,  the  analysis  used  fewer  elements  and 
was  less  costly  to  run  than  one  using  a  full  three-dimensional  model.  For  this  analysis,  the  two-dimensional  region  was 
subdivided  into  triangles  and  quadrilaterals.  In  an  axisymmetric  analysis,  nonaxisymmetric  loads  can  be  applied  to  an 
axisymmetric  structure  by  expressing  them  as  Fourier  expansions  in  the  azimuthal  angle.  However,  if  a  nonaxisymmetric 
loading  produces  nonlinear  effects,  it  may  introduce  non-physical  effects.  For  example  with  horizontal  accelerations, 
the  effect  of  the  interface  between  the  cylinder  and  displacer  is  nonlinear,  since  where  the  surface  is  in  compression, 
there  are  action  and  reaction  forces  between  the  cylinder  and  displacer,  but  where  the  surface  is  in  tension  there  are  no 
forces  between  them.  In  an  axisymmetric  analysis  the  interface  is  either  connected  or  disconnected  at  all  azimuthal 
angles.  It  was  assumed  that  these  effects  did  not  significantly  influence  the  results. 

The  cold  head  at  the  lower  end  of  the  first-stage  cylinder  fits  into  the  first-stage  cylinder  wall.  The  cold  head 
is  joined  to  the  cylinder  wall  by  a  weld  that  was  estimated  to  connect  the  lower  one-fourth  of  the  surface  between  the 
two  pieces.  See  Appendix  A  for  a  sketch  of  the  weld  detail.  For  horizontal  accelerations  it  was  assumed  that  the  cold 
head  fits  snugly  in  the  cylinder  and  the  relative  normal  motion  of  the  surfaces  was  constrained  where  they  were  in  contact 
but  not  welded.  For  the  region  where  the  surfaces  were  in  compression  this  was  a  physically  realistic  constraint,  but  for 
the  region  where  these  surfaces  tended  to  separate,  it  was  not.  The  finite  element  analysis  showed  that  for  vertical 
acceleration,  these  surfaces  tended  to  separate  around  the  cylinder,  so  no  constraint  on  relative  normal  motion  was 
applied  for  this  case. 
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A  motor  drives  the  displacers  through  a  scotch-yoke  mechanism.  To  avoid  modeling  the  driving  mechanism, 
the  analysis  was  limited  to  the  situation  in  which,  at  the  time  of  the  shock,  the  displacers  were  at  the  lowest  point  in  the 
cycle  so  that  each  was  resting  on  the  bottom  of  its  cylinder.  For  vertical  acceleration,  maximum  load  on  the  cylinders 
and  thus  maximum  stress  occurs  when  the  displacers  are  resting  on  the  bottom  of  the  cyhnders.  For  horizontal 
acceleration,  the  moment  of  inertia  about  the  attachment  of  the  refrigerator  to  the  vacuum  vessel  is  largest  when  the 
displacers  are  at  the  bottom  of  the  cylinders  and  also  produces  maximum  stress  in  the  cylinders.  It  was  also  assumed 
that  there  was  sufficient  slack  in  the  connection  joining  the  two  displacers  that  they  could  be  assumed  to  have  been  not 
connected. 

The  displacer  of  each  stage  in  this  design  of  the  Gifford-McMahon  refrigerator  contains  a  regenerator,  a  cavity 
filled  with  heat-absorbing  material  that  contributes  to  the  mass  of  the  displacer,  but  does  not  contribute  to  its  strength. 
Structural  mass  contributes  significantly  to  the  stiffness  of  the  object  and  is  represented  by  finite  elements;  when  this 
material  is  assigned  a  mass  density,  the  finite  element  program  automatically  computes  its  contribution  to  the  mass 
matrix.  Non-structural  mass  makes  at  most  an  insignificant  contribution  to  the  structural  stiffness  and  therefore  is  not 
modeled  with  finite  elements.  Non-structural  mass  is  combined  with  the  structural  mass  of  adjacent  elements  and  in  this 
analysis  was  included  by  an  appropriate  increase  in  the  density  assigned  to  the  material  of  these  elements.  The  structural 
part  of  each  displacer  was  modeled  with  finite  elements  leaving  a  cavity  for  the  regenerator.  The  mass  of  the  regenerator 
material  was  determined  and  a  layer  of  finite  elements  on  the  surface  of  the  cavity  was  assigned  a  density  that  accounted 
for  the  layer’s  mass  and  the  non-structural  mass  of  the  heat-absorbing  material.  For  horizontal  acceleration,  the 
non-structural  mass  was  assigned  to  the  layer  of  elements  on  the  side  of  the  cavity  to  represent  horizontal  inertia  and 
for  vertical  accelerations,  the  non-structural  mass  was  assigned  to  the  layer  of  elements  on  the  bottom  of  the  cavity  to 
represent  vertical  inertia. 

Normal  motion  on  the  surfaces  between  the  displacers  and  cylinders  was  constrained  and  tangential  motion 
between  these  surfaces  was  not  constrained.  With  axisymmetric  modeling,  this  assumption  implied  that  for  horizontal 
loading  the  displacer  pushed  on  the  cylinder  wall  if  its  relative  motion  was  toward  it,  but  it  also  implied  that  the  displacer 
pulled  on  the  wall  if  its  relative  motion  tended  to  be  away  fi'om  the  wall.  The  first  implication  was  physically  realistic; 
the  second  was  not.  For  vertical  loading  these  constraints  caused  no  problem. 

The  material  property  data  were  obtained  for  the  three  structural  materials  represented  in  the  analysis;  copper 
and  type  304  stainless  steel  which  form  the  cylinders,  and  a  cotton-phenolic  composite  which  forms  the  displacers. 
Where  temperature  dependent  data  were  available,  data  for  the  4K  operating  temperature  of  the  refrigerator  were 
chosen.  Table  1.  gives  material  property  values  for  three  temperatures  and  for  the  values  used  in  the  finite  element 
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Table  1.  Material  Property  Values 


Stainless  Steel  304 

Used 

19K 

30K 

300K 

Reference 

Youngs  Modulus  -  GPa 

199. 

199. 

200. 

192. 

Figure  B.4.P 

Shear  Modulus  -  GPa 

76.5 

76.5 

78.9 

73.6 

Figure  B.2.P 

Density  -  Kg/M^ 

7851. 

Copper 

Used 

4K 

30K 

300K 

Reference 

Youngs  Modulus  -  GPa 

139. 

139. 

139. 

128. 

Table  1.2^ 

Shear  Modulus  -  GPa 

51.7 

51.7 

51.7 

47.6 

Table  1.3^ 

Density  —  Kg/M^ 

8960.- 

Cotton-Phenolic  Composite 

Reference 

Youngs  Modulus  -  GPa 

6.27 

NEMA  Grade:  LE 

Shear  Modulus  —  GPa 

2.14 

General  Electric  Grade:  1841 

Density  -  Kg/M^ 

1335. 

Mil  Spec  15035FBE 

analyses.  Because  densities  for  stainless  steel  and  copper  were  used  to  determine  the  mass  of  the  finite  elements,  and 
the  dimensions  of  the  structure  are  the  dimensions  at  room  temperature,  the  densities  at  room  temperature  gave  correct 
mass  values.  The  elasticity  data  for  stainless  steel  and  copper  were  taken  from  the  Cryogenic  Materials  Data  Handbook .  ^ 
The  value  of  Young’s  modulus,  E  was  obtained  from  the  “Annealed”  curve  in  fig  B.4.i.  The  equation  G  =  E/2(l  +  v) , 
withv  =  0.3  determined  the  shear  modulus,  G.  The  data  for  the  composite  was  taken  from  General  Electric  data  sheets 
for  cotton-phenolic  composite  designated  NEMA  Grade:  LE,  General  Electric  Grade:  1841,  MilSpec  15035FBE.  From 
Young’s  modulus,  the  preceding  equation  determined  G  also  for  the  composite.  In  some  cases  for  the  three  materials 
values  and  significant  figures  were  estimated  from  plots. 


FINITE  ELEMENT  MODEL 

An  axisymmetric  structure  is  a  structure  that  is  rotationally  symmetric  about  an  axis  and  whose  shape  can  be 
generated  by  rotating  a  two-dimensional  plane  area  around  the  axis.  A  finite  element  model  for  an  axisymmetric 
structure  can  be  made  by  dividing  the  generating  plane  area  into  finite  subregions  or  elements.  The  boundaries  of  the 
finite  elements  consist  of  lines  or  curves  through  points,  called  nodes,  in  the  plane  area.  Thus  the  nodes  represent  grid 
circles  generated  by  the  points  as  they  are  rotated  about  the  axis  of  symmetry,  and  the  elements  represent  solid  rings 
generated  by  the  rotated  areas.  The  finite  elements  used  to  model  the  refi’igerator  are  the  NASTRAN  axisymmetric 
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triangle  TRIAAX  with  three  nodes,  and  the  axisymmetric  trapezoid  TRAPAX  with  four  comer  nodes."^  These  elements 
approximate  the  continuous  displacement  field  inside  the  element  by  linear  functions  that  interpolate  between  displace¬ 
ments  at  the  nodes.  Both  axisymmetric  and  non-axisymmetric  loads  can  be  applied  to  these  elements  to  produce 
axisymmetric  and  non-axisymmetric  displacements  and  stresses.  Non-axisymmetric  loads  and  displacements  are  ex¬ 
panded  in  harmonic  series  and  the  coefficients  in  the  series  become  the  finite  element  unknowns.  The  vertical  accelera¬ 
tion  was  applied  as  an  axisymmetric  load  and  the  horizontal  acceleration  load  was  applied  as  a  non-axisymmetric 
one-term  cosine  series.  For  axisymmetric  loading  or  at  most  a  horizontal  acceleration,  the  displacements  are  of  the  form: 


radial 

Ato 

“f  Ati  cos  (}) 

azimuthal 

Aai  sincf) 

axial 

Azo 

-1-  Azicostf) 

For  axisymmetric  loading,  the  terms  involving  ^  do  not  appear. 

Since  the  cylinder  walls  are  thin  compared  to  their  length,  a  large  number  of  elements  was  needed  for  the  model. 
To  represent  bending  adequately,  four  elements  through  the  thickness  of  the  cylinder  walls  were  used.  To  obtain  a 
well-conditioned  system  of  equations,  the  ratio  of  the  longest  to  shortest  side  should  be  less  that  1.5.  These  two 
requirements  determined  the  number  of  elements  in  the  cylinder  walls.  Constraints  linking  grid  points  on  the  cylinder 
walls  with  corresponding  points  on  the  displacers  also  determined  the  number  of  elements  in  the  displacers. 

A  refined  mesh  accommodated  the  larger  stresses  expected  in  the  region  of  the  weld  that  joins  the  two  stages. 
A  mesh  generation  program  written  for  this  model  produced  the  major  part  of  the  finite  element  grid.  The  refined  mesh 
in  the  region  of  the  weld  and  the  elements  representing  the  curved  surface  at  the  top  of  the  second  stage  was  added  by 
hand.  The  complete  model  contained  23350  elements  connecting  25011  grid  points.  The  figures  in  Appendix  B  show 
the  finite  element  grid  used. 


FINITE  ELEMENT  ANALYSIS 

This  section  provides  a  brief  outline  of  the  finite  element  method  as  it  applies  to  this  problem.  The  description 
follows  the  NASTRAN  Theoretical  Manual,  Section  5.11^^  which  can  be  consulted  for  the  details  of  the  method.  The 
finite  element  analysis  approximates  the  continuous  displacement  field  in  a  deformed  structure  by  a  set  of  displacement 
values  at  a  finite  number  of  locations  or  nodes  in  the  structure.  It  determines  the  set  of  these  displacements  that  produces 
the  minimum  potential  energy  of  the  system  consisting  of  the  deformed  structure  and  the  forces  producing  the  deforma¬ 
tion.  The  displacement  values  at  the  nodes  of  the  finite  element  model  are  the  components  of  a  displacement  vector  {u}. 
The  approximation  produces  a  system  of  linear  equations  which  is  solved  for  the  displacements  at  the  nodes. 
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The  energy  associated  with  a  displaced  structure  is  the  sum  of  the  strain  energy  U  of  the  elastic  deformation 
and  the  displacement  energy  W  of  the  loads. 

U  +  W  = 

where  in  the  volume  S  occupied  by  the  structure 
o  is  the  stress  field, 

e  is  the  strain  field, 

u  is  the  displacement  field,  and 

F_4  is  the  external  force  field. 

The  preceding  equation  must  be  expressed  as  the  finite  element  matrix  equation  [K]  {u }  =  —  {F^ } ,  a  system  of  linear 
equations  where  [K]  is  a  matrix  that  represents  the  stiffiiess  of  the  structure,  { }  is  the  vector  of  acceleration  forces, 
and  {u  }  is  the  vector  of  the  unknown  nodal  displacements.  The  analysis  was  made  using  the  NASTRAN  finite  element 
program. 

The  system  of  linear  equations  that  determines  the  displacements  at  the  nodes  of  the  finite  element  model  is 
the  sum  of  systems  of  equations  for  the  individual  elements.  The  following  description  of  the  analysis  is  given  for  one 
triangular  element  with  its  three  nodes  designated  a,  b,  and  c.  The  derivation  for  the  trapezoidal  element  is  similar,  with 
four  nodes  in  place  of  the  three.  In  the  analysis  of  axisymmetric  structures  with  nonsymmetric  loads,  the  displacements 
and  loads  are  represented  as  Fourier  expansions  in  series  of  sine  and  cosine  terms  of  the  azimuthal  angle  (J).  The  degrees 
of  freedom  are  harmonic  coefficients  of  the  displacement  components.  The  special  case  of  acceleration  loading  needs 
only  the  zeroth  and  first  terms.  With  a  vertical  axis  of  symmetry  the  zeroth  harmonic  displacements  for  vertical 
acceleration  are: 

u  =  Uo(r,z)  radial 
V  =  0  azimuthal 
w  =  Wo(r,z)  axial 

and  the  first  harmonic  displacements  for  horizontal  acceleration  are: 
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u  =  Ui(r,  z)  cos  4)  radial 
V  =  Vj(r,  z)  sincf)  azimuthal 
w  =  Wi(r,  z)cos(t)  axial 

The  linear  functions  of  the  radial  and  axial  coordinates  r  and  z  that  approximate  the  displacement  field  in  each  element 
are 


u  =  pi  +  Pjr  +  P3Z 
V  =  P4  +  Psr  +  PeZ 
w  =  Pt  +  PgF  +  P9Z 


For  a  triangular  element,  the  components  of  displacement  at  its  three  nodes  a,  b,  and  c  are:  Ua,  Vg,  Wa,  Ub,  Vb, 
Wb,  Uc,  Vc,  Wc  and  these  displacements  are  included  among  the  components  of  the  vector  {u}  and  are  the  components 
of  a  vector  designated  {u*}. 


The  vector  of  strains  for  the  n  =  0, 1  harmonics  are: 


5Un 

dr 

dz 

r  Vn  T  r 

dr 

-  ^^u  + 

J*  ' 


d\ 


dr 


-  -w 
dz  r  " 


“  [H^(r,z)](P} 


The  matrix  [H^(r,  z)]  is  obtained  using  the  Unear  approximation  of  the  displacement  field  and  is  a  function  of  the 
position  coordinates  r  and  z  and  the  harmonic  parameter  n.  {p}  is  a  vector  with  components  Pi, . .  ..Pg.  Now  the  vector 
{p}  must  be  expressed  in  terms  of  the  displacement  vector  {u}. 

Since  the  individual  harmonics  are  independent,  each  can  be  analyzed  independently.  Subsequently,  the 
subscript  n  designating  the  displacement  harmonic  is  suppressed  in  favor  of  a  subscript  i  designating  the  displacement 
at  the  ith  node  of  an  element.  The  linear  approximation  applied  to  the  three  nodal  points  (ri,Zi)  for  i  =  a,  b,  c  of  the 
triangular  finite  element  produces  the  transformation: 
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“0 

fPl  +  PjFj  +  PjZj' 

Vi  ■  = 

-j  p4  -f  PjTi  +  P^Zj 

[Wij 

IPt  +  Ps^i  + 

Using  these  transformations,  the  vector  {p}  and  the  vector  {u®}  with  components  Ua,  Va,  Wa, . . .  ,W(;  we  define  the  inverse 
of  the  matrix  [Hpu]>  determined  only  by  the  coordinate  locations  of  the  element  nodes,  by; 

{u‘}  =  [Hp-‘]{P}  so  that  (P)  =  [HpJ{u‘} 


The  strain  field  is  now  approximated  in  terms  of  the  vector  of  displacements  by: 

{«„}  =  [H"p(r,z)][HpJ{u‘} 


and  the  stress  vector  is  given  by: 


{o}  =  [E]{£} 

where  the  matrix  [E]  depends  only  on  Young’s  modulus  and  the  Poisson’s  ratio.  Now  in  each  element,  the  stress  field 
is  approximated  by: 

{o'}  =  [E][H."p(r,z)][HpJ{u'} 

where  {o*}  is  the  restriction  of  the  stress  vector  to  the  elements  nodal  points. 


Since  the  components  of  the  vector  {u*}  are  are  included  among  the  components  of  {u},  the  transformation 
matrix  [HpJ  can  be  extended  to  a  transformation  of  {u}  and  the  expressions  for  {e*}  and  {o®}  can  be  extended  to  {e} 
and  {a}.  In  terms  of  nodal  displacements,  the  strain  energy  of  the  displaced  structure  is  approximated  by: 

U  =  i  ^  {u}^  f[HpJ^H^(r.z)l^[E][H^(r,z)][Hpu]dV  {u} 

elements  J 

=  i{ur[K]{u} 

The  stiffness  matrix  [K]  and  the  displacement  vector  {u}  form  the  first  term  of  the  finite  element  matrix  equation.  The 
force  vector  remains  to  be  obtained. 


The  mass  of  an  element  having  area  A,  radial  distance  to  its  centroid  ro,  and  density  q  is  M  =  2jtQroA.  Each 
node  of  a  triangular  element  is  assigned  one-third  of  the  element’s  mass  and  each  node  is  assigned  the  sum  of  the  masses 
for  all  elements  attached  to  the  node.  If  the  mass  assigned  to  a  node  is  mj  where  the  index  j  ranges  over  all  nodes  in  the 


9 


finite  element  model  and  if  the  node  is  subjected  to  an  acceleration  A,  the  component  of  force  applied  to  the  node  is 
mj^.  The  potential  energy  of  acceleration  forces  on  the  displaced  structure  is 

W  =  {F^}  •  {u} 

Where  {F^}  is  the  vector  of  acceleration  forces  on  the  structural  nodes.  For  equilibrium,  the  derivatives  of  the  total 
potential  energy  with  respect  to  the  individual  displacements  must  vanish 

^(U  + W)  =  [K]{u}  +  {F^}  =  0 

This  is  the  system  of  linear  equations  that  forms  the  finite  element  matrix  equation  and  is  solved  for  the  displacements 
at  the  nodes  of  the  finite  element  model.  In  most  cases,  this  ^stem  requires  adjustments  to  reflect  boundary  conditions, 
constraints,  etc.^ 


ANALYSIS  RESULTS 

Three  types  of  analyses  were  made  using  the  finite  element  model.  Those  of  the  first  type  were  analyses  of 
applications  of  constant  lOOg  accelerations  applied  separately  in  the  vertical  and  horizontal  directions.  Displacements 
and  stresses  produced  by  these  loads  were  plotted.  The  second  type  was  an  analysis  of  an  application  of  a  time-depen¬ 
dent  vertical  acceleration  at  the  top  of  the  first  stage.  For  the  time-dependent  acceleration,  the  maximum  stress 
components  were  identified  and  the  accelerations  at  the  lower  end  of  the  first  and  second  stages  were  plotted.  The  third 
type  was  a  determination  of  natural  vibration  frequencies  of  the  refrigerator  with  the  top  of  the  first  stage  mounted  on 
fixed  flange.  The  refrigerator  was  oriented  with  its  axis  of  symmetry  vertical  and  with  the  first  stage  at  the  top  (see 
Figure  3.  in  Appendix  A).  The  acceleration  loads  were  applied  at  the  top  of  the  first  stage  where  the  refrigerator  is 
attached  to  the  vacuum  container.  The  goal  of  this  work  was  to  determine  whether  stresses  in  the  cylinders  of  the  two 
stages  would  remain  below  the  elastic  limit,  so  only  stresses  in  the  cylinders  were  computed  and  examined. 

Constant  Acceleration  Case 

In  the  NASTRAN  model  for  the  constant  acceleration  case,  the  acceleration  load  was  applied  by  constraining 
the  top  of  the  first  stage  and  applying  a  constant  gravitational  acceleration  of  lOOg  to  all  elements  in  the  model.  Two 
loading  cases  of  lOOg  acceleration  each  were  applied  to  the  model.  The  first  load  was  applied  in  the  vertical  direction 
and  the  second  in  the  horizontal  direction.  The  resulting  displacements  and  stresses  were  plotted  for  each  case. 
Appendix  B  contains  examples  of  these  plots.  Table  2.  shows  the  maximum  absolute  value  of  the  stresses  for  each 
loading  case  and  each  stress  orientation.  The  stresses  observed  were  less  than  the  yield  stress,  785  MPa  for  cold  rolled 
tempered  304  stainless  steel.^ 


10 


Table  2.  Maximum  Time-Independent  Stresses 


Vertical  Acceleration  -  100  g 

Radial  Stress  -  MPa 

11.4 

Axial  Stress  -  MPa 

36.7 

Circumferential  Stress  —  MPa 

16.5 

Shear  Stress  -  RZ  —  MPa 

11.1 

Horizontal  Acceleration  —  100  g 

Radial  Stress  —  MPa 

37.8 

Axial  Stress  —  MPa 

143.5 

Circumferential  Stress  —  MPa 

83.1 

Shear  Stress  -  RZ  —  MPa 

36.8 

Displacements  for  grid  points  on  the  cylinders  were  plotted  and  printed.  Table  3.  gives  the  maximum  displace¬ 
ments  for  the  two  loading  cases.  The  maximum  radial  displacement  for  the  horizontal  case  occurs  at  points  on  the  bottom 
of  the  second  stage  lying  on  the  plane  defined  by  the  center  of  the  cylinder  and  the  displacement  vector;  the  maximum 
azimuthal  displacement  occurs  on  the  cylinder  90  degrees  fi’om  these  points.  The  values  of  these  maximum  horizontal 
displacements  are  equal. 


Table  3.  Maximum  Time-Independent  Displacements 


Vertical  Acceleration 

-lOOg 

Component 

Value  -  mm 

Location 

Radial 

.0001 

Bottom  of  Stage  1 

Axial 

.0172 

Bottom  of  Stage  2 

Horizontal  Acceleration  -  100  g 

Component 

Value  —  mm 

Location 

Radial 

.5715 

Bottom  of  Stage  2 

Azimuthal 

.5715 

Bottom  of  Stage  2 

Axial 

.0787 

Bottom  of  Stage  2 
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Table  4.  Natural  Frequencies  of  the  Refrigerator 


Frequencies  (Hz) 

Period  (sec) 

1324 

0.00076 

2400 

0.00042 

3108 

0.00032 

Natural  Frequencies  of  the  Refrigerator 

An  eigenvalue  solution  of  the  finite  element  matrix  equation  gives  the  natural  frequencies  of  the  modeled 
structure.  Table  4.  shows  the  natural  frequencies  of  the  refrigerator  model  and  the  corresponding  periods  of  the 
vibrations. 

Time-Dependent  Acceleration  Case 

The  time-dependent  acceleration  data  were  taken  from  the  curve  “Lightweight  Test,  Top  Blow”  in  Figure  15 
in  Bradley  et  al^ .  This  curve  gives  the  vertical  acceleration  at  the  test  machine.  The  maximum  amplitude  of  this 
acceleration  was  103  g.  In  the  NASTRAN  model  for  this  case,  the  time-dependent  acceleration  load  was  applied  to  the 
top  of  the  first  stage.  Since  transient  analyses  produces  a  large  amount  of  data  for  each  node  selected,  stresses  were 
selected  to  be  computed  only  at  nodes  in  regions  that  showed  high  stress  values  in  the  static  case.  Table  5.  gives  the  time, 
location,  and  value  of  the  maximum  stresses  at  these  nodes. 

The  time-history  of  the  accelerations  at  the  bottom  flanges  of  the  two  stages  were  computed.  Figure  2.  shows 
the  vertical  acceleration  applied  at  the  flange  of  the  vacuum  vessel  and  the  vertical  components  of  the  computed 
acceleration  at  the  two  stages.  The  radial  and  axial  components  of  the  computed  accelerations  were  small  compared 
to  the  vertical  components.  The  curves  in  Figure  2.  show  that  the  input  acceleration  changes  slowly  relative  to  the  period 
of  the  lowest  natural  frequency  and  therefore  the  response  follows  the  applied  acceleration  closely.  The  frequency 
content  of  the  applied  acceleration  does  not  approach  resonance  and  the  maximum  response  amplitude  was  smaller  than 
that  of  the  constant  acceleration. 


CONCLUSIONS 

The  maximum  stress  produced  by  the  analysis  for  vertical  steady  state  shock  loading  was  36,7  MPa  (from 
Table  2,).  This  value  was  5  percent  of  yield  for  tempered  stainless  steel.  The  maximum  stress  produced  for  horizontal 
shock  loading  was  143.5  MPa  which  is  19  percent  of  yield.  This  stress  was  closer  to  the  yield  value,  however  it  is  believed 
that  any  horizontal  shocks  experienced  by  this  equipment  will  be  much  less  than  the  assumed  lOOg.  The  maximum  stress 
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found  in  the  time-dependent  analysis  was  18.4  MPa  (from  Table  3.)  which  is  2  percent  of  yield.  Comparison  of  the 
natural  frequencies  in  Table  4.  with  the  curves  in  Figure  2.  shows  that  the  time  varying  acceleration  applied  to  the 
structure  changed  slowly  compared  to  the  period  of  its  natural  frequency.  Therefore  the  time  varying  acceleration  did 
not  excite  a  resonance  in  the  structure  and  produced  smaller  stresses  than  the  constant  acceleration  of  similar  maximum 
amplitude.  As  a  result  of  these  analyses,  given  the  anticipated  values  of  horizontal  and  vertical  shock  during  testing,  it 
is  expected  that  the  refrigerator  will  not  experience  damaging  strains. 
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APPENDIX  A  -SKETCHES  OF  CYLINDERS  AND  DISPLACERS 


Figure  3.  is  a  sketch  of  the  cross  section  of  the  first-  and  second-stage  cylinders;  it  gives  the  dimensions  used  in 
the  finite  element  model.  Figure  4.  is  a  sketch  of  the  displacers  showing  the  dimensions  used.  The  detail  of  the  welded 
joint  between  the  two  stages  is  shown  in  Figure  5.  The  depth  of  the  weld  is  assumed  to  equal  1.65  mm.  The  surfaces 
that  are  welded  are  assumed  to  be  in  contact  in  the  region  that  is  not  welded.  Because  of  expected  stress  concentrations 
near  the  top  of  the  weld,  the  finite  element  model  was  refined  in  the  region  of  the  weld. 
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Note:  dimensions  in  mm 


22.22 


76.20 


(a)  First  Stage 


13.97 


13.97 


29.92 


25.40 


(b)  Second  Stage 


Figure  4.  First  and  Second  Stage  Displacers  (approximately  to  scale) 
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Figure  5.  Detail  of  the  Welded  Joint 
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APPENDIX  B  -  PLOTS  OF  DISPLACEMENTS  AND  STRESS  LEVELS 


Plots  were  made  showing  the  displacements  and  stresses  produced  by  the  finite  element  analyses.  This  appendix 
includes  plots  of  the  undeformed  and  deformed  finite  element  meshes  and  color  stress  level  plots  for  the  displacement 
components  having  the  largest  stresses  for  each  of  the  shock  loading  cases  that  were  run.  Figure  6,  shows  the  unde¬ 
formed  finite  element  model  for  the  cylinders  for  the  two  stages  and  the  deformed  plots  for  the  horizontal  and  vertical 
shock  loading.  Individual  elements  are  not  visible  in  this  figure  because  of  the  large  number  of  elements.  Since  the 
maximum  displacements  in  the  horizontal  and  vertical  cases  were  of  the  order  of  0.02  mm  and  0.57  mm  respectively,  the 
plotted  displacements  are  greatly  exaggerated  to  be  visible.  Figure  7,  shows  enlargements  of  the  displacement  plot  for 
the  vertical  acceleration  for  the  region  in  which  the  two  stages  are  joined,  the  lower  end  of  stage  one  and  the  upper  end 
of  stage  two.  These  enlarged  plots  show  the  refinement  in  the  finite  element  grid  for  the  region  near  the  weld.  Even 
though  the  displacements  were  very  small,  it  is  evident  that  the  largest  distortions  were  in  the  vicinity  of  the  weld  and 
the  connection  between  the  first  and  second  stages. 

Color  plots  showing  the  stress  levels  in  the  cylinder  walls  were  made  for  all  of  the  stresses  shown  in  Table  2. 
For  each  component,  one  plot  shows  the  stress  levels  in  the  complete  finite  element  model.  The  other  plot  is  an 
enlargement  of  the  region  which  connects  the  two  stages.  This  region  contains  the  welded  joint  that  joins  the  two  stages 
and  the  curved  section  that  fairs  the  lower  end  of  the  first  stage  into  the  wall  of  the  second  stage.  These  two  locations 
produce  some  of  the  higher  stress  values  experienced  in  the  analyses.  Figure  8.  and  Figure  9.  are  plots  of  the  axial  stress, 
the  largest  stress  components  for  the  horizontal  and  for  vertical  shock  accelerations.  For  the  horizontal  acceleration. 
Figure  8.  shows  that  the  stresses  in  the  vertical  direction  were  largest  where  the  bending  moments  are  largest,  at  the  top 
of  the  first-  and  second-stage  cylinders.  These  axial  stresses  were  due  to  compression  on  one  side  of  the  cylinder  and 
tension  on  the  other.  There  were  also  axial  stress  concentrations  near  the  top  of  the  weld  that  joins  the  first  and  second 
stages.  Figure  9.  shows  axial  stress  levels  for  the  vertical  acceleration  case.  In  this  axisymmetric  case,  axial  stress  was 
due  to  stretching  or  in  some  locations  compression  which  was  uniform  around  the  cylinder.  The  largest  stress  concentra¬ 
tions  were  at  the  top  of  the  weld;  positive  on  one  side  and  negative  on  the  other. 
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Undefonned 


Deformed  by  horizontal  shock  load 


Figure  6.  Undeformed  Plot  and  Deformed  Plots  for  the  Two  Shock  Loads 
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Figure  7.  Undeformed  and  Deformed  Displacement  Plots  for  Lower  End  of  Stage  One 


20 


Axial  Stress  -  Ibs/in^ 
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Figure  8.  Stress  Plot  for  Horizontal  Acceleration 


Axial  Stress  -  Ibs/in^ 
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Figure  9.  Stress  Plot  for  Vertical  Acceleration 
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